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Abstract 

We develop the functional renormalization group formalism for a tensorial group field theory with closure 
constraint, in the case of a just renormalizable model over (7(1)®®, with quartic interactions. The method allows 
us to obtain a closed but non-autonomous system of differential equations which describe the renormalization 
group flow of the couplings beyond perturbation theory. The explicit dependence of the beta functions on the 
running scale is due to the existence of an external scale in the model, the radius of ~ (7(1). We study 
the occurrence of fixed points and their critical properties in two different approximate regimes, corresponding 
to the deep UV and deep IR. Besides confirming the asymptotic freedom of the model, we find also a non¬ 
trivial fixed point, with one relevant direction. Our results are qualitatively similar to those found previously 
for a rank-3 model without closure constraint, and it is thus tempting to speculate that the presence of a 
Wilson-Fisher-like fixed point is a general feature of asymptotically free tensorial group field theories. 
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1 Introduction 


Among the many results obtained with the renormalization group, two of the most famous are certainly the 
asymptotic freedom of QCD and the non-trivial IR behavior of scalar field theories below the critical dimension, 
due to the presence of a Wilson-Fisher fixed point. However, the two phenomena are hard to find in the same 
system: the non-trivial IR behavior of QCD is due to confinement, not to an IR fixed point, and the marginal 
couplings of scalar field theories do not go to zero in the UV, running instead into Landau poles. One example 
of coexistence is provided by gauge theories with Nf massless fermions, with Nf close to the critical number 
above which asymptotic freedom is lost: in such case one finds the so-called Banks-Zaks fixed point [T]. However, 
such a fixed point is not a generic feature of gauge theories, it exists only for a restricted range of parameters. 
More recently, a new class of field theories has emerged which appears to enjoy asymptotic freedom in quite some 
generality, and for which we will show here another case of its coexistence with a Wilson-Fisher fixed point. Such 
theories are known as Tensorial Group Field Theories (TGFTs). 

TGFTs are a particular case of Group Field Theories (GFTs) |5], a class of field theories defined on a group 
manifold, whose non-local interactions give a complex cellular structure to the Feynman’s graphs of the pertur¬ 
bation theory. Such a structure is particularly interesting for quantum gravity, because it provides a tool for 
summing over manifolds. In their simplest form, and historically their original one, they appear as generalizations 
of matrix models [3] to higher dimensions, in the form of tensor models |1] . The most fruitful form of these models 
are the colored tensor models, for which an expansion in 1/N, similar to the one for matrices, has been constructed 
by Gurau and his collaborators 0 [SJ [3 HI IS] • Such results boosted a fast expansion in this direction (e.g. [TU] 
references therein). Proper GFT models can be seen as enrichments of tensor models with group-theoretic data. 
They originated in the context of topological field theory m, and they were further developed because of their 
links with loop quantum gravity miini and spin-foam models M- Lastly, TGFTs can be seen either as an enrich¬ 
ment of the colored tensor models by the addition of group variables, or alternatively as GFTs whose interactions 
show the same structure and the same unitary invariance as the interactions of the colored tensor models. Such 
particular models will be studied in this paper, by means of a specific example. 

From the quantum gravity point of view, one of the main challenges of GFTs is to understand how the quantum 
degrees of freedom organize to form a geometric structure which can be identified with a semi-classical space-time 
(a challenge that GFT share with loop quantum gravity m)- In this quest, the physics of phase transitions 
seems to be a promising route. In the most widely admitted scenario m, and according to recent works nzms], 
these phase transitions should correspond to the condensation of quanta corresponding to the fields involved in 
the GFTs. This phenomenon, very similar to the Bose-Einstein condensation, implies the spontaneous acquisition 
of a nonzero value by the mean field (something similar has been explored recently also in tensor models [191 [20]). 
In this approach, the help of the renormalization machinery, which perfectly matches with the GFTs formalism, is 
invaluable m- The renormalization of GFTs |22| and TGFTs |23112411251126j has been studied for several years, 
and the TGFTs (with non-trivial propagator) are the only models found so far which allow to build rigorously a 
renormalization program, a success mainly due to the loosening of the ultra-locality of usual GFTs, and to the 
tensorial nature of the interactions. This leads to a new notion of locality, called “traciality”, which enables one 
to clearly define the contraction notion of a high graph [23]. For now, several approaches have been followed, 
especially perturbative ones, and some characteristics of these models start to emerge. One of them is their 
asymptotic freedom, which has been discovered by Ben Geloun and his collaborators by perturbative calculations 
|23| ■ and which seems to be a generic characteristic of TGFTs j^. Another aspect is the presence of non-trivial 
fixed points, which started being explored only recently by means of the Functional Renormalization Group (FRG) 
|291130] and the e-expansion m- 

The FRG is a widely used formulation of the Wilsonian RG, with a vast range of applications |321 1531 [3311331 
IMIHZIISS]. For many purposes its most convenient formulation is in terms of the effective average action, which 
satisfies what is often called the Wetterich equation [311110]. The main advantage of the FRG is that it allows 
to explore different approximation methods that do not necessarily rely on having a small coupling constant, and 
thus are particularly interesting in the search of non-trivial fixed points. It is precisely for such reason that we will 
employ here the FRG machinery, following the path opened by its application to matrix models (Hj and TGFTs 

[20J|3g. 

In this paper, we will further develop the FRG approach to TGFTs, by studying the case of a just renormalizable 
Abelian TGFT [^. After briefly reviewing TGFTs in section i and the FRG construction at the beginning of 
section we will obtain, by the end of section a closed system of differential equations which describe the 
renormalization group flow of the model. In sections and we will study the occurrence of fixed points and 
phase transitions in two different approximate regimes, corresponding to the deep UV and deep IR, or to the limit 
in which the radius of the group manifold goes to infinite and zero, respectively. We will summarize and discuss 
our findings in section leaving the connection to the perturbative calculations for the appendices. 
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2 Tensorial Group Field Theories 


2.1 Generalities 

A tensorial group field theory on d copies of the compact group G is built along the lines of uncolored tensor models 
of rank d jH]- The theory is defined by a choice of Gaussian measure dfici'^pj V') of interaction Sint['<pj V'] the 
generating function 




( 1 ) 


where C refers to the covariance of the free theory. The fields '0 and ip, as well as the sources J and J are functions 
on G‘^ that take values in C, 

ip:G'^^C J : ^ C, 


and the notation (•, •) means 


{J,ip) := J d‘^[g]J{gi,...,gd)ip{gi,.-,gd), (2) 

where d[g] stands for the Haar measure on the group G being considered. Obviously, if G has dimension D, the 
TGFT is a field theory in dimension D x d. We denote by T and T the set of fields ip and ip, respectively. The 
true peculiarity of TGFT models lies in the non-local form of the interactions contained in Smt, which are built as 
follows. The convolution of the fields follows the outline of index contractions from the interactions of the uncolored 
tensor models jH]: group elements are contracted (i.e. identified and integrated over the group) pairwise, such that 
an element appearing as the i-th argument of a '0 is contracted with an element appearing as the i-th argument 
of a 0. A pattern of contracted fields is sometimes referred to as a bubble. What makes the bubbles non-local is 
both the restriction to pairwise identification and the fact that a given field can contract different arguments to 
different fields. While the interactions are non-local on the group, the pairwise contraction allows the definition 
of a new locality principle, known as traciality, essential for the renormalization. In the lattice gauge field theory 
language, traciality means essentially that, for all the leading order graphs, fiat holonomies imply necessarily 
trivial connection [21]. In the case that G is the abelian group U{1), this amounts to imposing invariance under 
a formal unitary symmetry, the limit for N going to infinity of the uncolored tensor models’ invariance. 

The interactions are formed by all the possible index contractions compatible with the invariance, i.e. 

they are formed by the pairwise contraction described above, with no insertion of external functions or operators 
on the group. We refer to these interactions as invariant traces'. 

'•= ^ A6Tr5[V;,'0], (3) 

b 

where b labels each contraction possible, involving an arbitrary number of fields, and is associated to the coupling 
Af,. Those contractions are graphically represented by bipartite colored graphs, whose vertices, black or white, 
are associated to the fields ip and ip, respectively, and whose lines correspond to the contracted field indices. For 
example, the graph of figure below, 
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Figure 1: Example of invariant bubble in dimension 6 
corresponds to the following rank-6 interaction: 


52 , 53 , 54 , 55 , 56 ) 0 ( 5052 , 53 , 54 , , 55 , 56 ) 0(5050537 505056 ) 0 ( 51 , 5050505056 )- (4) 


Other examples of bubbles are given on figure]^ for various ranks. There is obviously only one bubble constructible 
out of two fields, represented in the last graph of figure This is the only local interaction in the usual sense, 
it corresponds to a mass term, and it is often included directly in the covariance. All the other interactions arise 
from the multiplication of fields at different points in G®'^. 
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Figure 2: Examples of interaction bubbles in dimension 5, 4 and 6 

In the TGFT models, a specific symmetry is often imposed at the level of the field variables, known as the 
closure constraint, which reads 


tp{9i,---,9d) =ij{gih,...,gdh) ^heG. (5) 

This gauge symmetry is part of the physical enrichments added to purely tensorial models in the GFT context 
and it can be interpreted as geometrical data. Its origin can be found in the spin-foam foundations and LQG 
requirements of group field theories HIH], and recent works have shown that it implies some nice properties in a 
purely field theoretical context. In particular, it has been pointed out that it might favor asymptotic freedom, see 

m- 


The choice of the covariance C must respect this invariance and must enforce it at the level of the Feynman graphs. 
The usual choice is the following, given here in Schwinger (or parametric) representation. 


f' p p-\-oo ^ 

JG Jo 


where m > 0 is the mass parameter of the theory and is the heat kernel on G, verifying 

d 

— Ffa = J^gKa, 

Ag being the Laplacian operator on G. 


( 6 ) 

(7) 


As in the standard quantum field theories, the calculation of perturbations to the Gaussian model can be 
evaluated thanks to the Wick theorem, and any fV-points Schwinger function, Sn, can be evaluated as a sum over 
Feynman graphs. These graphs, of which figure ULis an example, are formed by all possible Wick contractions, 
represented by dashed lines, between all couples 



Figure 3: A typical Feynman graph in d = 4 with four external lines 
A 7V-points Schwinger function can thus be written as 

( 8 ) 

where s(0) is the dimension of the automorphism group of the Q graph (the symmetry factor), b label the vertices 
of Q and Ag is the graph amplitude. The expression of the amplitude can be obtained by using the elementary 
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properties of the heat kernel. We find: 



In this formula, s{e),t{e) are the source and terminus of the oriented edge e, Ce/ designates the incidence matrix, 
being 0, +1 or -1 given that the line e belongs or not to the boundary of /, the sign depending on the relative 
orientation, and C(Q),J^(Q) and J^ext{Q) are respectively the sets of lines, internal and external faces. Note that 
in the language of lattice gauge theory, the integrals over the group variables h^, in the amplitude ([^ can be 
understood as a sum over discrete connections on the cellular complex defined by the Feynman graph Q. From 
this point of view, the product IleGa/holonomy around the face /, and the symmetry of 

the Feynman amplitudes (where |I^(t/)| is the number of vertices in Q) coming from the invariance under the 
transformations 

he — t ks[e)hek^f^^-^ , ky (z G, (10) 

is a discrete version of gauge invariance. 

In this paper we study the case G = U{1) and d = 6, that is, a six-dimensional Abelian model, with interactions 
limited to the degree 4 melonic bubbles of the same form as that of figure All the interactions of that kind, for 
each possible choice of the solitary line color, are included in the model with the same coupling constant. In this 
model, the heat kernel has a simple expression in momentum representation 

iF„(6») (11) 


and the integration on the group d[h] in the propagator’s definition shows a so that the expression 

of our model’s propagator is written, in momentum representation 


C{v,p’) 


= 5 : 


p,p '; 




9 ’ 


( 12 ) 


where S{x) := stands for a Kronecker delta. The interaction part can be written as: 


[[de}j (13) 

where bi designates the bubble of figurej^whose solitary line is of color i and (8)T^ —)■ C is a super-tensor 

of order four, invariant under unitary transformations]^ As mentioned before, these interactions are chosen to be 
melonic, meaning that their Gurau degree |7] vanishes: vj = 0. In momentum space we have: 



^Fi!p2.P3.P4 ~ ^Pli,Pii^P2i,P3, 


n •P2j ^P3j ,Pij ■ 


(14) 


As explained in |S], the Gurau degree characterizes the colored graphs associated to the interaction bubble, and 
interestingly, it plays a role similar to the genus for the matrix models in the 1/N expansion of the colored tensor 
models. Melons appear as the leading order graphs in this framework, and it can be proved that they are dual 
to topological spheres. Their names comes from the leading order graphs of the colored tensorial models, called 
melons in j^. Our model is uncolored, in the sense of [3], but our Feynman graphs admit a natural colored 
extension, and one can show that their corresponding Gurau degree vanishes for the melonic graphs. 

Because of the non-local structure of the interactions, the renormalization becomes quite non trivial. However, 
recent works have proved that a renormalization scheme is possible, thanks to the traciality property of the 
divergent graphs (see |24| 1. coming from the tensorial structure of the interactions. The power counting of our 
model was carried out using multi-scale analysis, and it has been proven | 23 ] that the divergent degree oj{Q) for 
an amplitude Ag is given by: 

coiQ) =- 2 L{g) + F{g) - R{g), ( 15 ) 

^That is, for a fixed size Ai of a tensor, Wb) g inv({/(fV)®'^) 
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where L{Q) and F{Q) designate respectively the number of lines and faces of the graph mentioned, and R{G) the 
rank of the incidence matrix defined below equation (|^. It can be checked that our model is just-renormalizable, 
with divergent graphs of degree 2 and 0 requiring only mass, field and coupling renormalization. More interestingly, 
it can be proved that all the divergent graphs are melonic, in the sense explained before. Hence, for our model, 
the melonic sector contains all the graphs that need to be renormalized. 


2.2 Canonical dimension of a tensorial bubble 

In order to investigate the influence of the dimension d on the occurrence of fixed points in the FRG analysis, 
we need to work with dimensionless quantities, as in usual quantum or statistical field theories. In the present 
case, there exists an external scale L (the radius of ~ C^(l)) which we have set to 1, thus implicitly making 
all quantities dimensionless in the usual sense. In fact, as pointed out in [2^, and as usual in the presence of 
an external scale |42| . this is the reason why we will obtain non-autonomous RG flow equations. However, what 
matters for our purposes is how such quantities scale with the cutoff, which we refer to as scaling dimension. This 
dimensional notion appears quite naturally in the perturbative calculations. In Schwinger regularization, which 
introduces a cutoff 1/A^ in the lower end of the integral over a appearing in ([^, we get, at the dominant order in 
A and in dimension d: 


Siioop.oo =AA‘^-4Xi 

r[1Uoo=24(-A + A^A‘^-6iF2), 


(16) 

(17) 


where Ki and K 2 are two numerical constants independent of A, A or m (see appendix]^. The exponent of A in the 
above expressions is in general a universal number, meaning that it does not depend on the chosen renormalization 
(for example, we find the same exponent by regularizing the theory by means of a cutoff on the momentum). We 
therefore define the (canonicaQ scaling dimension [X] of a quantity X in such a way that by redefining X = XA^^l 
we obtain a homogeneous expression in A for the renormalization of X. The equation 0 then shows that the 
scaling dimension of A is [A] = 6 — d, so that both summands have the same dimension. By using this result in 
(16), we find [Si loop, 00 ] = 2, implying that [to] = 1, since Ei loop, 00 gives the radiative corrections of the mass term. 


We can easily check the coherence of these definitions at all orders, noting that the exponent of A appearing 


in (161 and ( [17| ) is nothing else but the divergence degree uj of the corresponding graph. The divergent degree is 
given by 0, and it can easily checked |23] that the most divergent graphs, the so-called melons, verify: 


F{g) - R{g) = {d- 2){L{g) - v(g) + 1), 


(18) 


where V{g) is the number of vertices in the G graph. The equation (181, together with the combinatorial relation 
HG) = 2V{G) - No^t{g)/2, leads to 


uj{g) = {d-6)V{G)+ {d-2) + 


4-d 


No^tiG) 


(19) 


We limit ourselves to the interactions that we can find in the initial theory. A IPI graph with four external legs 
has, a priori the same dimension as the coupling constant A. That is to say: 


(d - 6 + [A]) (R(e) - 1) + (d - 6) + (d - 2) + X 4 


4-d 


= 0 , 


implying [A] = 6 — d. The same argument applied to an IPI function with two external legs justifies [to] = 1: 

= 2. (20) 


(d-6 + [A])R(0)+ (d-2) + i^x2 


3 FRG approach for TGFT 

In this section we apply the Functional Renormalization Group formalism to the TGFTs introduced above. First 
we derive the Wetterich equation [Ml SO] in this context, and later we use it to extract the flow equations in a 
simple approximation. For general reviews and applications of the FRG we refer to jSH ESI El ESI ISa EZl ESI • 

^At a non-trivial fixed point the scaling dimension will in general be anomalous, i.e. different from the canonical one. 
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3.1 Effective average action for TGFTs 

Starting from the generating functional defined by the equation Q, we define the following deformation: 

where we have added to the action an IR cutoff or “momentum-dependent mass term” AS'g, chosen ’’ultralocal” in 
the momentum representation, and formally defined as: 

A5,[^,:= (^,= Y, Rs{p)Tp% (22) 

peZ'^ 


where is the Fourier transform (or Peter-Weyl decomposition) of the field: V’(^) = X)peZ‘i ^■ As usual, in 
order to make sense of the partition function we assume that a UV regulator is also present, e.g. a sharp cutoff on 
the momenta |p| < A (for the usual norm: \p\ = y/p ■ p). In practice, one often works in the limiting case A —> oo, 
because the Wetterich equation is well defined in that limit (although not its path integral origin). 

The cutoff function Rs{p) is a positive definite function chosen so that: 

• Rs{p) > 0 for all p G and s S (—oo, -foo). 


lims_,._oo Rsil^ = 0> implying: 


Zs=-o,[J,J]=Z[J, J]. 


(23) 


This condition ensures that the original model is in the family (211. Physically, it means that the original model 
is recovered when all the fluctuations are integrated out. 


• lim^-^inA Rs{p) = + 00 , ensuring that all the fluctuations are frozen when e® = A. As a consequence, the bare 
action will be represented by the initial condition for the flow at s = In A. 


• For — oo < s < In A, the cutoff Rg is chosen so that 

Rg{\p\ > e®) « 1, (24) 

a condition ensuring that the UV modes \p\ > e® are almost unaffected by the additional cutoff term, while 
i?s(|p| < e®) ~ 1, or Rs(|p| < e®) » 1, will guarantee that the IR modes \p\ < e® are decoupled. 

• ^Rs{p) < 0, for all p & 1/ and s G ( —oo,+oo), which means that high modes should not be suppressed 
more than low modes. 


As it stands, (211 defines an infinite-dimensional deformation of the original partition function. However, the role 
of the precise cutoff function, chosen to satisfy the above requirements, is secondary with respect to the role of 
the parametric dependence on s. From a Wilsonian point of view, the former corresponds to a choice of coarse 
graining scheme, while the latter corresponds to the coarse graining scale. We are primarily interested on the scale 
dependence of the theory, and therefore we take the point of view that a specific cutoff function has been chosen, 
and view (21) as a one-parameter family of theories|^ We thus obtain a one-parameter family of free energies. 


Wg :=logZg[J,J], 


(25) 


and by their Legendre transform, a one-parameter family of effective actions, collectively called effective average 
action. To be more precise, the effective average action F^ is defined as: 

+ {^,Rgcf) = {JA) + {^,J) - Wg[J,J], (26) 


where the source J is to be expressed as a function of the effective mean field (p via the solution of 


SJ 


(27) 


The previous properties concerning the cutoff term ASg mean, at the effective average action level, that: 


^In principle, physical quantities, such as critical exponents, are independent of the coarse graining scheme (see for example the 
universality of the one-loop beta function for marginal couplings, which we discuss in appendix^jl, but approximations generally spoil 
this property. Scheme dependence is thus an important issue, which has been greatly developed into the art of optimization |43| . but 
we will not discuss it further here. 
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• = Sint, so that when all the fluctuations are frozen, the effective average action coincides with the 
initial “microscopic” action. 

• rs=_oo = r, meaning that when all the fluctuations are integrated out, the effective average action coincides 
with the full effective action. 


The definition (21) of the modified action can be usually interpreted as a modification of the propagator. For a 


standard Gaussian measure, with usual kinetic action, this property is obvious, but it is not so obvious for our 
model, for which the covariance has been defined via equation ([^ rather than by an explicit kinetic action]^ To 
prove this point, we use the Wick-theorem. Consider the following covariance: 


Cs(p,p') = 




Pt 


(28) 


and the two Gaussian integrations: 

CMp) ■■= J dficAT,f)f{p)T{f) 

Js{p,p') := / 


(29) 

(30) 


where d^c is the normalized Gaussian integration: f dfic = 1- In a first step, we will prove that (x Cs- This 
comes obviously from the Wick theorem. Using the derivative representation of a Gaussian integral: 


Js{p,p') = exp 


T,T=0 


(31) 


ax yOip dtp 


d 

dx 
d 

~ Ex^ 


T,T,x=0 


-Trln(l-|-C(fls-|-a;L)) 


= Cs{p,p') X det 
a;=0 L O 


where the elements of the (super-) matrix L are defined as Nsxt, consider the following 

Gaussian integral: 


N 

J(:= J dptcAT,f)l[fAPj)TAp')- 


(32) 


i=i 

Using the Wick theorem, and the previous result, we obtain: 


= EH (/ dpc'{T,f)f^«^^\p^^(^))TAp'S 

■Kn j ^ 

= En(d®t ^ |d/rc(r,T)e-<W> xT-««(p^„(,))T^(p;)Y 


TTiV J 


where is the permutation group of iV elements. Now, using the Wick theorem, the big parenthesis can be 
written as: 


/, := det 


C_ 

c: 


^ n! 


dpc{T, (p,„(,) )T^ (pj) 

dpc{T,f) X {f,RsTAf^-(^Ap..U))THPj) 


c 

Cs 


= det 


EE 


_1 '\n-p+p 


(- 1 ) 


nl 


n\ p\{n — p)\ 

t' 

X ((T,i?,r)p)^ ((T,i?,r)"-pr-«(^)(0;^(,.))r^(0;) 


(33) 


C,NV ' 


®Alternatively, we can define our covariance C in its momentum representation, directly via equation ( |12[ l. In this case, we see that 
because of the projector the definition of the inverse C*”^, which usually appears in the kinetic action, is slightly subtle and 

we prefer to present a more formal derivation involving only C. 




















where the brackets {.)c mean the Gaussian contractions with covariance C, and {■)c,nv with subscript NV select 
only the “non-vacuum” contractions (i.e. involving the external fields only). The Cauchy decomposition formula 
allows to rewrite the double sums as: 

p \ y- ' C ri ’ ' 


The first term (^e~ ^ is nothing but the trivial Gaussian integral, i.e. the normalization of the Gaussian 

measure d^Cs'- 

^ g-Trln(l+CA) ^ ^ ^ ^ ^ 

implying finally: 


Reporting this result into the equation (331 for J' leads to: 


J's = Y.]l^s=det 


Tl’iV 3 




C,N.V 


En^ 

T^N j 


(34) 


Finally, because of the fact that the factor (e ^ is nothing but the vacuum contribution of the measure 

we find: 


N 


dfic. (T, f) Y[ y {p, )P (p') = det 
i=i 


N 


i=i 


(35) 


Not surprisingly, this result is exactly what is expected for a standard Gaussian integration with well defined 
kinetic action. Hence, it follows from the properties of the Gaussian integration that the definition (211, with 
modification of the action, can be interpreted as a modification of the covariance, i.e. 




(36) 


In fact, the flow equations for (211 and (361 differ for vacuum terms due to the det(Gs) implicit in the measure of 
(36) (see (35)), but as usual such vacuum terms are completely unimportant. 

Because of the gauge invariance of our model, it is useful to define the projector into the gauge invariant field 
P as: 


P-.r-^G (37) 

P P[ilj]{0i,...,ed) = J + v^ ...,94 + 11), 

and the invariant field subspace as G = ker(P — I) C T. Interestingly, because of the gauge constraint in the 
propagator, the gauge constraint is implemented, at the graph level, at each colored line hooked to a black or white 
vertex. Hence, the mean field (p, defined by ( |27[ ) lives in G. The gauge symmetry is then dynamically implemented 
by the propagator, and contaminates each A^-points functions : Gjv € Gi^iG. 


For the 1-point function (p, this property can be proved as follow. From the general expression ([^, cp can be 
expanded as: 

E ^(n(-^^))-^e( 0 ), ( 38 ) 

Q connected ^ \bGQ / 

where Ag{9) is a Feynman amplitude of order V{Q) with one external point. A typical graph contributing to this 
expansion is depicted on the Figure below, and its explicit expression has the following structure: 


Ag{9)= / n4 n 


(39) 


j=l 


l '=2 
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where for convenience we have used of the propagator Cg defined by (28 1 unlike of the propagator C. From its 
definition, the propagator lives in G 0 G, because it verifies : 

PCsP = Cs. 

Hence, it follows that the three vector indices 02 ,^ 3,04 of the (super-)tensor - - are projected into the 

gauge invariant space, and this tensor is equivalent to : 


yv 


a)_ 


27r ^1,^2,^35^4 


(40) 


where 9i := {du + Pi, 02i + Pi, ...,6di + Pi)- The consequences are clearer in Fourier components: 


yv 


( 1 )_ 


/ 11 Ott Pl^P2,P3,. 


27T 

i—2 Pi 


gipi-Si X 


4 ^ , 5 


j=2 


(41) 


=nnd Eb* ■ 

Pi j=l 1=2 \k=l / 

The index structure of the bubble vertex implies that pik = P 2 fc V/c 7 ^ l,P 2 i = P 31 and pn = p^. But, because of 
the Kronecker delta, P 4 i = p^i, implying pn = p 2 i and finally pi = P 2 - Hence p— 1 satisfies the closure constraint: 
Pik = 0 , and the last term of the previous equation can then be rewritten as: 


W 


( 1 )_ 

)^ 2 1^3 5^4 




i=i 




2 = 1 


27r 01,82,83,64 


(42) 


meaning that P[p] = p, or that p G Q, thus proving our claim. 



Using the Legendre transform (26) and the definition (21), we obtain formally the so-called “Wetterich equation”, 
describing the evolution of the effective average action. We have: 


Proposition 1 (Wetterich equation) 

For a given cutoff Rg, the effective average action Fg satisfies the following partial differential equation: 


SgFg = ^ dsRsip) ■ [ri^) + Rg] \f^,p)6 



=Tr 


P- 


dgRg 


42 ) 


R.i 


(43) 


10 












where Tr is the ’’super-trace”, meaning the trace over the block-indices, and The presence of P restricts 

the sums over the subspace oflP which we denote by V and define by : V = {p G = 0}. 

Proof: The proof is rather standard Easn], but it is useful to review it because of the TGFT context and the 
appearance of the projector P in the equation. Applying the operator dg on the two members of the equation 
(26), we have 


dsTs = {dgj,fi) + {4>,dgj) - dgWg - {dsJ, - {fi,dsRs^). 


SJ 


6J 


(44) 


The term dgWa can be easily computed from the definitions (21) and (25): 


dgWg = - ^ dgRg{p) I 


5‘^Wg SWg SWs 


5 Jp 6 Jp 


SJp SJp 


(45) 


Using (27), many terms cancel in (44), and we are left with 


a^Ts = ^ dsRg{p) 


Deriving (27) with respect to J we find 


pGZ'i 

6^Wg _ 54> 
SJSJ ~ SJ’ 


SJpSJp 


while deriving (26) with respect to and then gives 

S'^T 


Therefore we obtain, in matrix notation: 


SfiScj) Sfi * 


(46) 


(47) 


(48) 


SJSJ 


S^Ts 

SfiSfi 


+ Rs 


= I 


(49) 


The second functional derivative is the 2-points function which, as the effective mean field (f = Ws^\ is 

gauge invariant: € G 0 G. Hence, the previous equation (49) means that T*-^^ -|- Rg is the inverse of 

on the subspace V only. Finally, using (49) in (46), and the fact that, because WU) is a gauge invariant matrix : 
€ G G G —)■ PW^^^P = we find the Wetterich equation (43). 


□ 


3.2 The truncation procedure 

One of the many advantages of the FRG is the possibility to make an approximation directly at the level of the 
effective average action. The latter is in principle a complicated functional, containing infinitely many terms, and 
therefore an exact solution of the Wetterich equation is generally beyond reach. One of the most used approxima¬ 
tion methods, besides perturbation theory, consists in truncating the space of functionals in which Fg is defined, 
retaining only a finite-dimensional subspace, or a simple infinite-dimensional one. This allows to obtain a closed 
and tractable system of differential equations, which can be analytically studied in some cases, or numerically in¬ 
tegrated for most of the cases. As an organizational principle for such truncations, one usually starts off by Taylor 
expanding the effective average action in powers of p. Truncating such expansion to some given order is expected 
to be a valid approximation in a context where we are interested in the long distance physics. Subsequently, one 
can proceed in expanding each order in powers of the fields, and truncate that expansion too, thus leading to a 
finite-dimensional subspace of functionals. The main idea is that by systematically expanding such subspace one 
should observe some convergence pattern in the values of physical quantities, such as the critical exponents, and 
thus justify a posteriori the approximation procedure. A neat example of this method is provided by the Ising 
universality class, associated to the Wilson-Fisher fixed point in three dimensions |44| . However, when applying 
the FRG in a new context, before reaching such advanced level of systematization, it is often necessary to start 
from the simplest possible approximation in such class of truncations, to work out the formalism and understand 
possible outcomes. In such case, the most usual choice essentially boils down to writing Fg exactly in the same 
form as the bare action we have in mind, but with coupling, mass and wave function normalization depending on 
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s. The latter is the approximation we will adopt here. 


More precisely, we will henceforth consider the following simple truncation: 


i [(1^, J dOfjiiO) [ - Zs + m^] (j){9) 


(50) 


where is the Laplacian operator on Ag and nis are the effective coupling and mass at the scale s, and 

Zs is the field strength normalization. 


Evaluating the second derivative of (501, we get, in momentum representation: 




ZsP + ml 


Pp,p' T 2As ^ p^SyTDyVp^^j;2,P3,PA^P3^pii 

{pi} 


. . . - is given in dl^, and we have introduced the notation W = V. , 

PI 1 P 25 P 35 P 4 ® 1_ 1 ’ ’ 


where wl* 


and 


Let us also define: 


and 


SymW^j^^2,P3:P4 •— '^pi ,P2 ,P3 ,P4, T ^P3,V 


^pi ,P2 ~ ^pi ,P2 ^ I 'y ^ 


Pi ■ 


PsiPiP ) ■ — 2 'y ) Pp,pi Pp' ,P2 Sym Wpi ,P2,P3,P4 '^P3 '^P' 


{pi} 


Ks{p) ■■= [Zsp'^ + ml + Rs 
The next step is the choice of the regulator Rs, for which we adopt Litim’s cutoff 

i?g(p) = Zg(e2*-p2)0(e2^_p'2)^ 


(51) 

(52) 

(53) 

(54) 

(55) 

(56) 


where 0 stands for the Heaviside step function. The special dependence on p and the inclusion of a wave function 
renormalization, will lead to substantial simplifications. Applying the derivative operator on i?s, we find: 


dsRsip) = |9gZg(e^^ -P^) + 2Zge^®|0(e^® -P^)- 


(57) 


Using the equations (511 and (571 in the flow equation ( |43[ ), we will obtain the renormalization group flow equations, 
which are the subject of the next section. 


3.3 Renormalization group flow equations 


Using our ansatz (501 in the Wetterich equation (431, expanding its r.h.s. 
order (jp, we find: 


of in power of p, and truncating at 


r.h.s. 


=E 

p&V 


dsRs 

Ks{p} 


Fs{p,p) , j ,2 Fs{p,p')Fs{p\p) 

" KPp) ^ ^ KMKsip') 


0{<p) 


The RG flow equations for coupling, mass and wave function renormalization can be obtained from this expansion, 
identifying the corresponding powers of cfxl)- 


• Equations for ms and Zs. Identifying the terms of order (jxp leads to: 


{dsml + dsZsP^)Tj^Tj; = 

pev 


- 2 As 

pGV 


{OsZsie"^^ -p^)+2Zse‘^^) 
(ZgB^® + mlP 


91.92 


(58) 
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We start the analysis with some considerations about the r.h.s.. Because of the form of the interaction matrices 
and Tp e G: 


p&'P qiA p,p'ev 


Defining p± to be the set (pi,... ... ,pd), and 


/(Pi) = I 


equation (591 can be rewritten as: 


E E = E f(p")Wp- 

P^'Pqi,q2 p^V 


In the same way, because of the closure constraint, 


E E = E f^P")TpTp. 

P^V qi,q2 p^V 


We can therefore rewrite (581 as 


E (dsml + dsZsP^)Tpfp = -2A« E 

p^V 


p£V \i—l 


f{pl) + df{p^)]TpTp 


(59) 


(60) 


(61) 


(62) 


(63) 


where 


and we have defined 


f{p^)={A + Bp^) e{e^^-p^), 
(9,Z, + 2 Z,)e 2 ^ 


A := 


B :=- 


(ZgC^® + m2)2 ’ 

dsZ, 


(ZgC^'^ + m ^)2 


We want to satisfy (631 up to order meaning that we want to expand 

d 

J2f{p^)+df{p^) = d(^m + f{0)) +p^ (/'(o)+d/'(o)) +o(r 


equate the coefficients of and p^, 


(64) 


(65) 

( 66 ) 


(67) 


= - 2 A 3 d (/( 0 ) + /( 0 )) , ( 68 ) 

= -2A, (r(0)+dr(0)) (69) 

and discard the rest. The expansion of /(p^) is trivial, as the step function is equal to one at all orders of the 
Taylor expansion, and we thus have /(O) = A and /'(O) = B. The expansion of f{pl) is instead slightly more 
involved. We can first write f{p^) as 

f[pl)=ASM) + BSM), (70) 

where 


Siik"^) :=^dik + Ep; “ k"^ -P±‘^), 




S2{k^) ■■= E(^^ +Pi)'^f^ + EP')^(®^* “ k^ 

p±_ \ l^i / 


(71) 

(72) 
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so that 


= -2A,d (A(1 + 5 i(0)) + B52(0)) 

dsZs = -2X, {AS[{Q) + B{d + S'^m)- 


(73) 

(74) 


Since we will be mostly interested in the large-s limit (the small-s case will be treated separately), we can ap¬ 
proximate the sums by integrals, replacing the Kronecker deltas by Dirac deltas. The support of the integrals is 
in the intersection of the hyperplane of equation k + = 0 and the (d — l)-ball of radius VNote 

that the Kronecker delta of the closure constraint can be rewritten as fc -|-pj_ • n = 0, where n = (1,1, ...1) is the 
vector with all components equals to 1 in Using the rotational invariance of our integral, we can chose one 

of our coordinate axis to be in the direction n. If we choose the axis 2 in this direction, our constraint writes as 
5{k -fp^ul) = 5{k -f p' 2 'ij d — 1) = S{k/\/ d — 1 + P 2 )IV d — 1, and we find the following integral approximation: 


S2{k^) 


1 




1 r dk"^ 


e 2 «- 


Vd^Ld- 1 


dk^ 
d-1 

d- 2 




=:Si{k^) 
dk"^ 


d-1 


I^d-2( e2«- 


dfc 2 

d-1 


=: ~S2{k^). 


(75) 

(76) 


where := 7 r‘^/^/r(d /2 -|- 1) is the volume of the unit d-ball. The integral approximations are easily expanded to 
yield: 


/i:=5i(0) = 


y/d^ 


IId- 2 e 


{d—2)s 


h := ^ 2 ( 0 ) = 

Vd - 1 d 

T _ c> /n\ _ ^ ~ o „(d-4)s 


h ■■= S^( 0 ) = - 


1 d(d - 4) 

Vd^ 2 (d- 1 ) 


^d-2& 


{d-2)s 


Therefore we obtain 


ds'm^g — — 2 As d (A (1 -|- Ii) -f B I 2 ), 
dgZg = —2As {AA B (d 14 )), 


which, using (65l-(66l and ([77|-(|^, translate into 


dgTni = -2\gZ, 


( 


2 ^d — 2 ^ds 


/d^ 


de 


-I- 2d(^ 


^d -2 ^ds I g2s ' 


(Z^e^® -I- TTlg)^ 


and 


where we have defined 


2d(d-2) AgUd_2e(®*-2)® 

(d - 1)^/^ (Zge2®m2)2 - 2Ag 

d 0 , ,c(d-2)s _ 1 _ rl 

5 

(rf_l)3/2 ^^d—2^ wU 


d 

Vs := ^log(Z,). 


(77) 

(78) 

(79) 

(80) 


(81) 

(82) 


(83) 


(84) 


(85) 


• Equation for Ag. In order to obtain the equation describing the running coupling constant flow, we must 
identify the correct melonic structure among the terms of order in the r.h.s. of the flow equation, i.e. in 

dgRs y2 Fs{p,p')Fs{p'jP) 

^^Kgip-) Kg^sm ■ 

The expression above involves a term of the form 

SymWj^ ,Pi ,P 2 .P 2 ’ 

giving three types of contributions, 

yyl - yvl ^ 

P,P ,Pl,P2' '^p,p' ,p[ ,^2 ’ 
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wt ^ ^ yvi ^ ^ 

Pl.P ,P,P2 p,p',p{,f2’ 

yvi ^ ^ yvi, ^ ^ 

Pl,P ,P,P2 Vl.p .P.^ 2 ’ 


all depicted in figure below. 


The connectivity of the original melonic interaction induces a selection rule for the interactions generated by 
the r.h.s of the flow equation. This selection rule is based on the study of the “ultra-local” version of the interac¬ 
tions depicted in figures [^, and|^. 


1 



Figure 5: Interactions involved in the r.h.s of the Wetterich equation. The dash-dotted line represents a contraction 
involving the operator K~^, and the dash-dotted line with a square represents a contraction involving the operator 
dsRs K-^ 

In order to define this selection rule, we need the help of the following definition, illustrated in figure 


Definition 1 (contraction operation). Let b an interaction bubble with L(b) dash-dotted lines. Let Lq = {h} C 
L{b) an ordered subset of dash-dotted lines in b. The contracted graph b/Lg is obtain from b by: 

• Deleting the line k S Lq and its two (black and white) end vertices and all the colored lines joining these two 
vertices. 

• Identifying the colored line linked to the deleted black vertex with the line of corresponding color linked to the 
white vertex. 

• Repeat for 



Figure 6: Contraction of a dotted line between two vertices. Some black and white vertices are labeled as Vi in 
order to facilitate the understanding of the picture. 

The purpose of this contraction procedure is just to allow us to obtain the connectivity structure of an interaction, 
and therefore we do not keep track of the black square (the insertion of the cutoff operator) in the definition above 
and in the figures. Of course we will keep track of that once we have identified the correct structures. As an 
example of contraction, consider the interaction depicted in [^, for which we obtain, after contraction over the 
two dash-dotted lines, the result depicted in figure]^ 
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Figure 7: Connecting structure of the interactions |^. 

We observe that the first contraction gives exactly a vertex of the original form, while the second contraction gives 
a disconnected vertex corresponding to the square of a mass-term, which is outside of our truncation, and thus 
will be discarded. As a second example, the case i = 1 and j = 2 for the last & case gives, after contraction of its 
two dash-dotted lines, the result depicted in figure showing that the interaction does not have the connectivity 
structure of a melonic bubble, and will also be discarded. 


1 



“ .i -> 

\ / 



1 



Figure 8: Connecting structure of the interactions with i ^ j 

Of all the possibles interactions, only those of type and & with i = j, and with i ^ j, respect the 
connectivity of the original melonic interaction. Therefore, we retain only such terms for the computation of dsXg, 
i.e. we make the approximation 

d 

,pi,P2^y^'^p',P,P3,P4 — ^ ^ 


_-k - wi. , 

P,P ,Pl,P2 P ,P,P3,P4 Pl,P ,P,P2 P3,P,P ,P4 


;(b 


;(b 


( 86 ) 


+ wyl, - - 

p,p'.P1.P2 P3,p,p' 


P4 


(Sym" 


W 




p,p ,Pi,P2,P3,P4 


after which we obtain the equation 


dsXg y^pi,p2,P3 


- - (9s^s(e^® — -I- 


p.p'gv 

^Sym^W 


-f nig)^ 


(87) 


7 ^-. 'T'-. 'T'-. 'T'-, 


' P,P',P1,P2,P3,P4 

As before, we want to satisfy this equation to the desired order in in this case meaning to zeroth order. The 
easiest way to accomplish that is to project the equation onto the field = Y[i=i leading to the equation: 


ddgXg = 4x1 E EX/ " 


P, |p|<e° p' 


(ZgC^® -I- m2)3 


Ep» j (Sym^W 


p,p',0,0,0,0 
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Using (141, and taking into account the closure constraint J2iPi = 0, we find - = nf=i 11^=1 


'^p,p',o,d = ni=i ^p,,pU ^ ^ consequence: 




m^W 


p,p',0,0,0,0 




Pi 


^p,p' • 


(88) 


i=l i=l 




As for the mass and anomalous dimension, the flow of As can be expressed in term of and S 2 , defined in 0 
and (721, giving 

+ e^^Si{0) - S2{0))dsZs+2e^^(^d + Sii0)^Z, 


5s As = 4A^ 


{ZgC^^ + TOg)^ 

or, in term of the integral approximations (l77| and ([?§: 


ff \ /I \2. 


5sZs -f 2 


Zg 


(Zse2® -1- m2)3 



(89) 


(90) 


3.4 RG equations for “dimensionless” parameters 

The RG flow equation obtained above describes the evolution of the couplings in our truncation with respect to 
the “scale-time” s. But, as in the quantum field theory or condensed matter models, a trivial contribution to the 
evolution comes from the “dimension” of the operator associated to these couplings. We define the dimensionless 
parameters Ag and frig as 

As = e*^®“‘^)®As ms = e^TOs, (91) 


using the results about the dimension of the “local operators” introduced in section 1.2. The flow equations for 
the dimensionless couplings are easily deduce from the previous ones. The equation for rjg is unchanged, and the 
equations for As and rris become: 


dsmj = —2ms ~ 2A, 


5sAs — {d — 6)As -f 4A( 


( 2 ^ + dgZs + Zs 

(Zs -k m2)2 

^ + I ^)5sZs + 2(de(2-rf)- + ^)Zg 

’ (Zs-km2)3 


Another more physical redefinition concerns the relations between the different renormalization parameters. 
The couplings As and mg which appear in the original model are not the effective mass and effective coupling at 
the scale s. Indeed, the evolution of the wave function renormalization Zg from 1 to an arbitrary value, modifies 
the effective couplings, which become : As = A”, fhg = \fZgfh'g. Hence, from the two previous equations, we 

finally deduce: 


5s m 


r2 

s 




^g(2-d)s^^s -k -k 

(1 -I- m”2)2 


(92) 


5sA: = (d-6-27?s)A”+4A: 


|'^g(2-d). ^ |^),7s + 2(de(2-'^)* + 


(1 


Tr2'l3 


(93) 


with: 


Vs = 


2d{d-2) 


^d-2^ 


(d- 1)3/2 (i^^r2)2_2Ar 


(d-l)3/2 


-k de(2-<i)« 


(94) 


Note that, similarly to what was found in j^l], the system of differential equations we just obtained is not au¬ 
tonomous, in the sense that even after switching to dimensionless couplings it maintains an explicit dependence 
on the RG scale s. This was interpreted in |29| (based on more general cases in which a similar phenomenon 
takes place, e.g. [1^) as consequence of the implicit scale set by the compactness of the group, i.e. the radius of 
~ C/(l), an interpretation corroborated by the recent results of |30| . where autonomous equations are obtained 
in the non-compact limit C/(l) —r M. 
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4 Large-s approximation and fixed points 

In this section we examine the flow equation in the large-s approximation, or in the UV limit. In this sector, the 
equations are much simpler, since they become autonomous, and thus it is possible to And flxed points, and their 
relevant or irrelevant perturbations. 


4.1 Large s approximation, vicinity of the Gaussian fixed point 


In the UV regime, corresponding to the large-s limit, the previous equations (92l, (93l, and (94l, for d > 2, reduce 
to: 


^r2 




-r2 r,\r 2r]s + 2d 




(95) 


= {d-6- 277,)A: + 4A 


r2 


^d-2 


dds + 2 


yjd-l (1 + ’ 


(96) 


2d{d-2)nd-2K 


(d-l)3/2(l+TO^2)2_2dr!^_2A^- ^ ’ 

These equation form an autonomous system, whose flxed points can be studied with standard methods. A trivial 
flxed point occurs for = A^ = 0, corresponding to the so called “Gaussian flxed point” (GFP). In order to 
study the stability of this flxed point, we expand the previous RGBs around this flxed point, at second order in 
the coupling constant. We obtain the following system (we limit the development of rjs to the first order in A^, 
because the anomalous dimension appears always as a quantum correction): 




d-2 


(d-l) 




y/d^ 


(d-iy 


(98) 


dsXl = (d-6)A”-4 


d{d-2) 


d-l 
2d{d - 2 ) 


- 2 


^d —2 \r 2 
^ ’ 


((f _ 1)3/2 

At d = 6 the coupling A” become marginal, and using O 4 = 7 r ^/2 we obtain: 

^r2 _ _n^r2 , -r2 _ il_2Tr' _ ^8 2Tr2 


(99) 

( 100 ) 

( 101 ) 




ds = 


5^5 

24^2 


( 102 ) 


(103) 


The minus sign in front of the r.h.s of equation ( 102[ ) means that the theory is asymptotically free. This result 
was found for several TGFT models iMiiMiissiiigr including our model |26[ [?7] , and it seems what it is a generic 
property of TGFTs models. Appendix give the one loop computation of the beta function, that we And in 
accordance with 102 and with the well known universality of the one-loop beta function (see Appendix [A| . 


For d > 6 , the theory is non-renormalizable, but the equations ( |98[ ) and (99) show the existence of a non-trivial 
flxed point for As and m at: 

- 1 -!)»/= 


A* = 


4 (d2 - 4d + 2)Od-2 ’ 


- *2 
m = 


d(d- 6 )(d-l) 


- 16d -f 16 


(105) 


d2-4d-f2 d(d-6)(3d-2)-4(d2-4d-k2)' 

For d = 6 the flxed point merges with the GFP, and for small e = d — 6 it is at a small value of the couplings. 


and therefore the use of the expansion (98 1 and (991 is justified (but note that the mass grows rapidly and 


diverges at d ~ 6.65). In order to analyze the stability of this flxed point, we compute the first derivative, at the 
point (m^*,A*) = ( 51 , 52 ) of functions Pi defined by (98) and (99), where the index i labels the couplings 
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A’") = (gi, 32 )- These derivatives build the “stability matrix” (3*^ = dj3i{g*)/dgj, whose expression for general 
d is not particularly enlightening. More important are its eigenvalues (the so-called “critical exponents”), which, 
at leading order in e, are —2 -|- 24e/7 and —e, respectively, meaning that we have two relevant perturbations (an 
operator is said to be relevant if it becomes important in the IR limit, irrelevant if it disappears in the IR limit, 
and marginal for vanishing critical exponent). Therefore, for small e > 0 we have asymptotic safety. Whether this 
survives at e = 1 is a question that we will not address here. 


For d = 6, (98) and (99) admit only the GFP, and the stability matrix reads: 


n* diag 
Pii 


_0 127r^ 

V5 


0 


0 


(106) 


with eigenvalues (—2,0) (for the GFP the critical exponents coincide with the canonical dimension of the associated 
operators), and eigenvectors, 

. 0 , 


h^ = 


ho = 



(107) 


The first direction is associated to a relevant operator (the mass), and the second one is marginal. To decide whether 
the marginal coupling is marginally relevant (i.e. asymptotically free) or marginally irrelevant (i.e. trivial) one 


has to go to next order in the expansion, as we have done in (102). 


Rather than expanding (98) and (99) in powers of the couplings, the flow can be studied numerically, and this 


is the subject of the next section (for the case d = 6). 


4.2 Non-Gaussian fixed point with d=6 


In d = 6, the autonomous system (95), (96), (97) has other non trivial fixed points than the trivial GFP. After a 


simple analysis of the equations, we And two fixed points, for the values: 


{Ai,mi} = 


r V5(43309 t 790141) -175 ± 0141 \ 


V 


21354847r2 


234 




(108) 


As for the Gaussian fixed point, we can study the stability of these fixed point, by computing the eigenvalues of 
the stability matrix, i.e. the critical exponents. Diagonalizing the stability matrix at the two fixed points, we And: 


A 


(+) diag 


4.86 0 

0 -0.9 


o(-)diag ^ f 262.8 0 

b I 0 7 


(109) 


( 110 ) 


showing that the Axed point {A* +, m* +} has one relevant and one irrelevant directions, whereas the other one has 
only irrelevant directions. 

Our result is very similar to the one obtained in [23] for a TGFT of rank-3 without closure constraint. As for 
that model, we And the existence of non-trivial Axed points (one of which carrying one relevant eigen-direction) 
in the critical dimension at which the model is asymptotically free. Note that, as we have seen in section |4.1[ 
moving inflnitesimally away from the critical dimension one loses asymptotic freedom, but a new non-trivial Axed 
point emerges from the Gaussian one. Such merging or splitting of Axed points near the critical dimension is a 


standard phenomenon, but it has nothing to do with the Axed points (108), which appear precisely at the critical 


dimension. Gontrary to that, this is not a standard phenomenon, and it has to be attributed to the non-local 
structure of the TGFT interaction. 

The numerical integration of the large-s flow leads to the phase diagram of flgure [9| w hich shows again a strong 
resemblance to the one found in |29| . Note that, because of our choice of truncation (|50| and parametrization, and 


in particular because of the choice of the regulator (56), the flow equations (95), (96), and (97) have singularities 
at A” = 5^/^(l-|-m''^)^/67r^. A similar singularity was found in |23|, where its origin is discussed in more detail. As 
a consequence of such singular line, the Axed point {A*_,m*_} is disconnected from {A*+,m*+} and the GFP. 

The non-Gaussian Axed point {A* + ,m*+} is reminiscent of the Wilson-Fisher Axed point. It appears as an 
IR Axed point associated to a broken phase, with positive coupling and negative mass. More precisely, the region 
above the irrelevant trajectory at the NGFP can be interpreted as a symmetric phase, with coupling and mass of 
the same sign, and the region below this irrelevant trajectory as a broken phase, with negative mass and positive 
coupling. 

Within the present truncation, our conclusion about the existence of such non-trivial Axed point is only valid 
in the large cutoA limit, i.e. in the UV limit. Therefore, the occurrence of a non-trivial behavior associated to the 
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broken phase, rests upon the survival of such fixed point in the IR, and this is why we will study the counterpart 
of the regime studied in this section, the small cutoff limit, in the next section. 



A. 


Figure 9: The flow diagram at large cutoff. The blue dot is the GFP, while the red and black are the non trivial 
ones, labelled with + and —, respectively. The black line corresponds to the singularity of the flow equations, 
because of which the shaded region is disconnected from the GFP. In green and brown are respectively the eigen- 
perturbation for the GFP and the NGFP. 


5 Small-s limit and comments 

When the cutoff becomes very small, e® < 1, all the sums obtained in |3.3| reduce to a single element, the zero 


mode p = 0. In this limit, for the sums defined in (711 and (721, we find S'! = 1 and S 2 = 0, independently of 


fc < e®. Therefore, from equations (|73|, (|74|, and (|89|), we find that the anomalous dimension vanishes, implying 

( 111 ) 
( 112 ) 


that we are free to choose Zg = 1, and that the beta functions becomes, in d = 6: 

dgml = -48Xg + 


dgXg = 56Xi 


(e^® -I- 


Using the rescaling mg = e®ms and As = e^®As leads to the following autonomous system: 

dgfhl = -2rhl - 48 

(1 -I- mj)^ 


dgXg — 4A,s 


56 


A 


(1 -I- TOg)^ ' 


In addition, as in the previous case, the system admits again a GFP, and a NGFP for the values: 

T* 49 12 

A = _ m* = -- 


(113) 

(114) 

(115) 


13718 19 

In order to study the stability of these fixed points, we compute the matrix /3ij as before. For the Gaussian fixed 
point, we find: 


— I ^ 

\ 0 -4 


(116) 


20 


























whose eigenvalues are (—4, —2), corresponding to two relevant eigen-directions. 

For the non-Gaussian fixed point, we find the matrix: 

/ M 17328 \ 

with eigenvalues (76/7,—2), giving one relevant and one irrelevant directions. Interestingly, the critical exponents 
are integers, a property found also in |29| . and which can traced back to the zero-dimensional nature of the theory 
in this regime (see discussion in |1S]). As for the large s limit case, the numerical integration of these flow equations 
give the Figure [T0| below. 

Note that in this IR regime the rescaling for mass and coupling does not correspond to the canonical one defined 
in section |2.2| Again, this is a consequence of the fact that at small-s our model is essentially a zero dimensional 
field theory, and our rescaling matches exactly with the one expected for a such theory, by fixing to zero the 
dimension of the (effective) action, and 1 the dimension of the mass parameter. Indeed, by direct inspection of 
our truncation ( |50[ ), we find [A] = 4. Because no phase transition is expected for such model, we have found an 
essential indication for the symmetric phase restoration in our model, as it should be on a compact spacej^ In the 
physics of phase transition, it is well known that phase transitions occur in the “thermodynamic” limit or, in other 
worlds, for a non-compact manifold. This limit, in our case, can be approached by restoring in our expressions 
the radius L of ~ U{1), here always fixed to 1, and by sending this radius to infinity. Such non-compact limit 
has been recently studied in m, confirming the fact that in such limit one is left with a picture identical to the 
large-s regime, and thus with a non-trivial phase transition. 

We have not extended our results in the intermediate cutoff regime, between small and large cutoff, but being 
the two asymptotic regimes qualitatively very similar to the ones in [29], we expect the intermediate regime to 
add no new qualitative insight with respect to the analysis performed there. 



A. 


Figure 10: The flow diagram at small cutoff. The color rules are the same as in figure Ordinary trajectories are 
in blue, eigen-perturbations of the GFP are in green, and eigen-perturbation of the NGFP are in brown. 


6 Conclusion 

In this paper we have developed the FRG formalism for a just-renormalizable model of TGFT on the compact 
group 17(1)®®, and we explored two regimes of the flow, for large and small cutoff, which are more suitable to 
analytical and numerical study. 

The lessons of our analysis can be summarized as follows. First, we have shown that the presence of the closure 
constraint can be incorporated in the FRG formalism without too much pain. We have thus obtained a Wetterich 
equation for this class of TGFTs. 

Second, from our analysis of a specific truncation, we have recovered in the UV limit the asymptotic freedom 
found perturbatively in ESI HZ] for this model, and more generally observed for several TGFTs models El Ha 

^For a related discussion on effective dimensional reduction and symmetry restoration see also m- 


21 









1151123]. Furthermore, in the same limit, we have found a non-Gaussian fixed point, of which we have studied the 
stability, and identified the critical directions. Interestingly, this non-trivial fixed point has qualitatively similar 
characteristics to the Wilson-Fisher one, and it is associated to a phase transition in the flow diagram. 

Lastly, in the IR sector, due to the compactness of the group manifold, we obtained an effectively zero¬ 
dimensional theory, and thus we have symmetry restoration in the deep IR. This is not a surprise, given the 
compact nature of the manifold. The compactness of space can be viewed as an additional regularization, to be 
removed in order to study the non-compact limit [Ml, or otherwise one can view the latter as the UV approximation, 
as we did here. In the non-compact limit one is left with just the large-s regime of our analysis (because s = In(fcL), 
if L is the radius of and k the dimensionful RG scale), and thus with a non-trivial IR fixed point associated to 
a phase transition between a broken and a symmetric phase. 

Overall, the above picture is strikingly similar to the one obtained in [^, and in part to the one in |30| . 
despite the models having essential differences. The model in is in three dimensions, whereas our is in six; 
the model in [23] has no closure constraint, whereas our does; and the model in [^ has a kinetic term linear 
in momentum, whereas our is quadratic. These are all important differences, and the fact that the final result 
is nevertheless so similar points towards the greater importance that lies within the main similarity between the 
two models: the non-local melonic interaction. We are therefore tempted to conjecture that such TGFTs with 
melonic interactions enjoy in quite some generality not only asymptotic freedom, but also a Wilson-Fisher-like 
fixed point. Of course, the approximation we employed here is the simplest one, and making our result more solid 
will require studying larger truncations. And more models will need to be studied in order to understand whether 
our conjecture applies to all asymptotically free TGFTs or only to some specific sub-class. Grasping also at a more 
qualitative level how such feature arises (e.g. similarly to how asymptotic freedom is understood to arise from the 
wave function renormalization taking place already at one loop |28| 1 is also an open question. We are just at the 
beginning of the FRG investigations of TGFTs and hopefully these and other aspects will be clarified very soon. 


Acknowledgment s 

We would like to express our gratitude to Joseph Ben Geloun, Daniele Oriti, and Vincent Rivasseau, for many 
useful remarks during the preparation of this work, and on the final draft. 


22 


A Universality of the one-loop beta function for A 

In this section we review a well-known argument explaining why the one-loop beta function does not depend on 
the choice of the regularization, and we give an explicit example of this result for the computation of the beta 
function with FRG method around the Gaussian fixed point using a different regulator than the one used in the 
main part of the paper. 

A.l A standard argument for one-loop universality 

Different computational scheme result in physical predictions being made in terms of different couplings. Suppose 
that in some scheme we have a dimensionless coupling A, whose flow equation reads 

d^X = biX^ + & 2 A^ + &3A^ + ... . (118) 

In a different scheme we have a different coupling A' with flow 

d,X' = b[X'^ + b'^X!^ + b'^X’'^ + ... . (119) 

At tree level the two couplings are the same, because tree level does not involve any ambiguity, thus the relation 
between the two couplings must be of the form 

X'= X + CxX^+ C2\^ + ... . (120) 


Now write 


ds\ = -^dtX 
oX 

= (1 + 2CiA + SCaA^ + .. .){biX^ + &2A^ + bsX^ + ...) 

= biX^ + { 2 Cibi + b 2 )\^ + (SCa^i + 2C162 + & 3 )A^ + ... 
= &iA '2 + 62 A '3 + (63 - 2C162 - Cfbi + C 26 i)A '4 + ... , 


where in the last row we used the inverse of ( 1201 , 


A = A' - CiA'^ - {C 2 - 2Cf)X'^ + ..., 


( 121 ) 


( 122 ) 


and we assumed that dtCi = 0. Therefore, under these assumptions, bi and 62 , i.e. the one- and two-loop 
coefficients of the beta function, are scheme independent. 

Note that if dsCi yf 0, the flow of A' is not autonomous, unless Ci = as -I- b, or unless the s-dependence of Ci 
is related to the beta function of another coupling. This is preciesely what happens in general within the FRG, 
which is a mass-dependent renormalization scheme |4ti| . However, when expanding the beta functions around zero 
mass, the above argument must hold true also for the FRG. 


A.2 Example with a different regulator 

As we mentioned in section the choice of a cutoff function corresponds to a choice of coarse graining scheme. 
Different choices thus affect the flow, but we expect that physical quantities, such as the critical exponents, will 
not depend on it. As we argued above, one such quantity is the one-loop coefficient of the beta function for the 
marginal coupling. To illustrate this point more explicitly, we choose here a different regulator: 


R^ip) 


ZsP‘ 


■ mz 


exp 


Z,p'^+ml 




- 1 


(123) 


and study the large cutoff limit (s ^ 1) of the flow equations around the Gaussian fixed point, in order to compare 
them with the ones obtained in section |4T| At the leading order inA^, i.e. discarding the derivative of Zg, the 
derivative of the regulator writes as: 


dsR'sip) 




Zgp"^ + Trig 


( 1 - 




Following the same procedure as in section [3.3| we And the three equations: 

ZflP^+17 


dginj. 


^ 4As 


^ Sym>Vp--o_g, 


(124) 


(125) 
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dsZ, ~ - 


4As d 


E 


^ ZsP^+Tl 

e-2'<e- 


'(E- 


Pi X SymW, 


p,p,k,k 


k^O 


dAs c. 44 E 


_ e 


.,o^(E4- 


ZsP + m^, 

For each of these equations the leading order can be easily extracted. Using the following sum formula : 


E^ 

n^Z 


, X—^ _ 1/3 + 277 

e 4a 


© E' 


n^Z 


/7r\i/2 

= U 


„ -iniiii , / 7rn/3 
1 + 2 ^ e “ cosh 


n>0 


the Fourier decomposition of and the distributional identity: 


we find: 


Using both the equations (1301 and (1311, we find for the effective coupling \s/Zl := A^: 

287r2 


dsK - - ^ 


a: 


( 126 ) 


(127) 


n 

(128) 

2 247r^ As 2s 

vs zf • 

(129) 

247r^ As 

(130) 

- 5©5 ^s' ’ 

47r2 A2 



(131) 


(132) 


in agreement with our results of section |4.1| and with the universality argument. We will see in the next appendix 
that a different computation of the one-loop beta function, using yet another regularization, gives the same result. 


B One loop computation of beta function in Schwinger regularization 

In this appendix we give some details about the perturbative calculation of the beta function at one-loop order. 

• 2-points function. The one particle irreducible (IPI) 2-points function divergences come from the melonic 
diagram depicted in figure El which contains five internal faces. From the definition of the divergence degree in 
section (2.2), we find that the divergence degree of such a diagram is equal to 2, implying that a wave function 


renormalization is needed to obtain a complete renormalized theory. Including symmetry factors, its expression 
writes as, in momentum representation: 




i=l J+i 


cp + vn? 


qi=Pi 


(133) 


In order to extract the divergent part of this sum, we use the Schwinger representation, which consist in a cutoff 
in the a integration. 




d/3 /■+°° 


*=1 




'l/A2 


Pi=qi 


Using again (127) and (128), we find, at leading order in A: 


/-A l^TT^A.n 2477^ A 2i 2477^ A , /,w2 

^teioip) = - 7^^ + —ln(A) + —— ln(A)p2 + C>(l/A). 


©5 


©5 


5V5 


(134) 


(135) 
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The first two terms correspond to a mass renormalization term, and we ignore them here. The third one is more 
interesting for us, as it gives the wave function renormalization Z \= \ Sz, with: 

5z=^^HA). (136) 



Figure 11: Melonic contribution to the IPI 2-points function at one loop order 


• 4-points function. The first one loop corrections to the IPI 4-point function involve three types of diagrams, 
but only those of the type depicted in figure contribute here. As for the tadpole of the previous calculation, 
they involve five internal faces, but a line more than tadpole graphs. Hence, the power counting of these graphs 
gives a logarithmic divergence. The structure of the IPI melonic 4-points function the following: 


p(4) melo 
^ Pl,P2,P3,P4 


^ [ - 4A + Ar(*)(p4.)] ^SymW«, 


(137) 


where: 


Ar«(p):=8A2 Y, 

Qj, jA* 


{<p + 



Making use once again of (128) and (127), we obtain 


ArWA(p)= l^ln(A) + 0(l/A) 


(138) 


(139) 


and, at the leading order in A: 


p(4) melo 
^ Pl,P2,P3,P4 


= -4A 


47r^A 


N — SymW-^. 

/ ^ 2 Pl,P2,P3,P4 



(140) 


Figure 12: Melonic contribution to the IPI 4-points function at one loop order 


•Beta function. The previous equations (140) and (1361 give the renormalization of the coupling and the 

(141) 


field at the one-loop order. The divergences are canceled by the rescaling: 


A = Z-^Z-^Xr 


, V' = Z^/'^'4)r 


with Z given at order A,, by (136), after replacing A with A^, and Z\ by: 

dTT^Aj. 


Z.=\- 


75 


ln(A) + 0(A7. 


(142) 
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Expanding the first relation (1411, we find, at order A^: 


A At’ 


47r^Ar 

V5 


ln(A) 


1 + 


487r=^Ar 

5^5 


ln(A) 


28 tt'^XI 


287r^Ar , ,. ’ 
1 + ln(A) 


5V5 


implying: 

Again this result is consistent with those obtained in section |4A] and in appendix [X| 


dX 287r^A^ 


(143) 


(144) 
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